Quantized coexisting electrons and holes in graphene measured using temperature 

dependent magneto-transport 
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We present temperature-dependent magneto-transport experiments around the charge neutrahty 
point in graphene and determine the amphtude of potential fluctuations s responsible for the for- 
mation of electron-hole puddles. The experimental value s « 20 meV is considerably larger than in 
conventional semiconductors which implies a strong localization of charge carriers observable up to 
room temperature. Surprisingly, in the quantized regime, the Hall resistivity overshoots the highest 
plateau values at high temperatures. We demonstrate by model calculations that such a peculiar 
behavior is expected in any system with coexisting electrons and holes when the energy spectrum 
is quantized and the carriers are partially localized. 
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Pristine graphene is a zero-gap semiconductor where 
the conduction and valence band touch at zero energy. 
This point is generally referred to as the charge neutrality 
point (CNP) with a zero charge carrier density. However, 
it has been shown experimentally, that for graphene on 
Si/Si02 substrates, electrons and holes coexist around 
the CNP, which is attributed to the presence of spatially 
inhomogeneous conducting electron- hole puddles. 

In this paper we investigate the coexistence of electrons 
and holes by means of temperature-dependent magneto- 
transport around the CNP. We present results in the clas- 
sical and quantum-Hall regimes. In the classical regime, 
the individual charge carrier concentration and conduc- 
tivity at the CNP increases with increasing temperature 
due to thermal activation but, as a consequence of poten- 
tial fluctuations, do not vanish even at the lowest temper- 
ature. From this temperature dependence, we determine 
the strength of the potential fluctuations responsible for 
the formation of electron-hole puddles, s « 20 meV for 
our samples. In the quantum Hall regime, the Hall resis- 
tance shows higher values than expected suggesting a de- 
activation of charge carriers with increasing temperature. 
Using model calculations, however, we demonstrate that 
this counterintuitive temperature dependence is straight- 
forwardly explained considering the quantized density of 
states of coexisting electrons and holes. 

The paper is organized as follows: In Sec.|TT]we shortly 
describe the samples an measurements. Sec. [Hi] is de- 
voted to the experiential magneto-transport in the clas- 
sical regime along with an interpretation of the results. 
In Sec. II VI we describe experimentally (subsection IIVI A,) 
and discuss theoretically (subsection II VB ) magneto- 
transport data in the quantum-Hall regime. The paper 



II. EXPERIMENTAL 

We have measured two field effect transistors made 
from single layer graphene (SLG) and bilayer graphene 
(BLG) deposited on Si/Si02 wafers and shaped into a 
1 /im wide Hall bar. Both flakes originate from nat- 
ural graphite and have mobilities ^ = 1 m~^ s~^ and 
jjL — 0.3 roT^ s^^, respectively. The total charge carrier 
concentration n is controlled by applying a gate voltage 
Vg on the conducting Si substratei^ 



(1) 



with a = 7.2 X lO^m^VV for a 300 nm thick Si02 gate 
insulator. The subscripts e and h indicate electrons and 
holes, respectively. 

Magneto-transport as a function of magnetic field, car- 
rier concentration and temperature were performed in a 
temperature range between 0.5 K and 250 K using a top- 
loading He-3 system and a variable-temperature insert in 
a 33 T Bitter magnet. 



III. CLASSICAL REGIME 

Fig. [IK a) shows Pxy{n) for the BLG sample at i? = 
0.8 T for different temperatures. Down to the lowest tem- 
peratures, the Hall resistivity does not exhibit quantized 
plateaus at this magnetic field, indicating a continuous 
energy spectrum. For all temperatures, the Hall resis- 
tance smoothly crosses zero at the CNP, which indicates 
that electrons and holes are present at all gate voltages<^ 
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FIG. 1: (color online) (a) Hall resistivity of the BLG sam- 
ple as a function of the total charge carrier concentration at 
B = 0.8 T for various temperatures, (b) Extracted carrier 
concentrations of electrons rie and holes for the same field 
and temperatures (see text for details). The solid lines show 
the calculated number of electrons and holes using the effec- 
tive density of states from Eq. Q and s = 19 meV. 



In this coexistence regime, the simultaneous contribu- 
tion of electrons and holes with equal mobility n to the 
Hall effect is given hyJ^ 



Pxy — 



B {n, - nh){l + {^lBf ) 
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Note that this expression originally derived for massive 
charge carries with mobility ^ = er/m* (r is the elastic 
scattering time and m* is the effective mass) is also valid 
for massless particles when defining the mobility more 
generally as the ratio between the average drift velocity 
of the carriers and the applied electric field. 

By solving the system of Eqs. ([T]) and ^ with re- 
spect to the variables rig and for each point of the 
experimental Pxy{Vg) dependence, we can determine the 
number of electrons Ue and holes rih independently;^ the 
results are shown in Fig. [TJb). 

At the CNP [n = 0) Ue and Uh are equal to a residual 
concentration uq, which is, according to Eq.[2l related to 
the slope of Pxy{n) at the CNP: 
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FIG. 2: (color online) (a) Electron and hole concentration no 
of the BLG sample at the CNP as a function of temperature. 
The solid line represents a fit to Eq. ((5)1 with s — 19 ± 2 meV. 
(b) Sketch of the random potential fluctuations in the sample 
yielding a broadened effective density of states shown in (c).— 



The temperature dependence of Uq is plotted in 
Fig. 121a). At zero temperature we find that a finite num- 
ber of both electrons and holes are present at the CNP. 
These electron- hole puddles originate from spatial poten- 
tial fluctuations within the sample represented schemati- 
cally in Fig.|2l^b) and lead to an effective density of states 
(DOS) around the CNP plotted in Fig.lHc)!^ 



D{E) = Do erfc(-£;/V2s)/2, 



(4) 



s is the amplitude of potential fluctuations. 

Though the question whether the puddles are caused 
by the impurity potential alone or by intrinsic ripples in 
graphene is still under discussion^ this phenomenologi- 
cal DOS is not sensitive to the deeper origin of the pud- 
dles. The calculated charge carrier concentrations using 
this effective DOS, shown by the solid lines in Fig. [TJb), 
indeed reproduce the experimentally observed behavior 
around the CNP. 

Directly at the CNP, the temperature dependence 
no{T) in a low temperature {kBT/s < 1) approximation 
is analytically found to bei^ 



no(T) = no(0) 



6 V s 



(5) 



and a flt of the data in Fig. |2{a) to Eq. |S] yields s = 
19 ± 2 meV. 

Analogously, the temperature dependence of the con- 
ductivity at the CNP can be described by:= 



CTxxiT) = aa:x{0) 
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(6) 



The fit of the experimental temperature dependence in 
the classical regime with Eq. (|6]) is shown in the inset 
Fig. 121 for the T trace. The fit yields s = 23 ± 2meV, 
which is close to the value s = 19 ± 2meV, determined 
from the fit to the Hall data. 
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FIG. 3: (color online) Zero-field conductivity of the BLG 
sample around the CNP for different temperature. The inset 
shows the temperature dependence of the minimum conduc- 
tivity at the CNP with the solid line as a fit to Eq. ((6l yielding 
s = 23 ± 2 meV. 



Remarkably, this strength of potential fluctuations is 
comparable to the thermal energy at room temperature 
and considerably larger than in conventional semiconduc- 
tor heterostructures.— Therefore, whereas quantum Hall 
localization in GaAs-based systems breaks down above 
30 K,— the strong potential fluctuations can induce a ro- 
bust localization of charge carriers in graphene when the 
Fermi energy is situated between two Landau levels, in 
particular around filling factor v = This allows the 
QHE to persist up to room-temperaturep^^*i^ 



IV. QUANTIZED REGIME 
A. Experimental Results 

Fig. mja) shows the measured Hall effect for BLG as 
a function of the total charge carrier concentration n at 
S = 30 T for various temperatures. For i? > 15 T the 
system is in the quantized regime and we observe plateaus 
at V = ±4, ±8 ... as expected for BLG. At the lowest 
temperatures, also a plateau develops a.t ly — S which 
makes the p^y curve slightly asymmetric. However, since 
the plateau is smoothed out for T > 74 K, we neglect 
this asymmetry in the following. 

At T = 0.5 K, the Hall resistance on the electron 
(hole) side reaches its maximum (minimum) value at the 
plateau and goes smoothly to the hole (electron) side as 
already observed in the non-quantized regime and this 
behavior can again be qualitatively explained by the two- 
carrier model. At high temperatures, the Hall resistivity 
overshoots the plateau v = ±4 before the beginning of 
the plateau and reaches values higher than /i/4e^ and 
Pxy{n) becomes steeper at the CNP, see Fig. EJa). We 
observe a similar behavior in the SLG sample, therefore. 



this effect has a generic nature and was also observed 
by others.^— Note that these overshoots in the Hall resis- 
tance only appears at high temperatures and is therefore 
not due to the (partial) splitting of the lowest Landau 
level as observed in e.g. Ref. 17i 

The temperature dependence of the slope of Pxy{n) at 
the CNP is explicitly shown in the inset to Fig. Ilja). It 
is opposite to the classical regime, which can be seen by 
comparison of the traces at _B = 30 T (quantized regime) 
and B = 0.8 T (classical regime), plotted together in the 
inset. 

According to Eq. ([3]), an increase in the slope with 
increasing temperature corresponds to a decrease of 
the charge carrier concentration, which appears to be 
counter-intuitive. However, the plateaus remain visible 
up to the highest temperature measured, i.e., the sys- 
tem remains quantized and the number of electrons and 
holes within one Landau level is fixed. In the quantized 
regime, the slope of the Hall resistivity depends not only 
on the number of charge carriers at the CNP, but also 
the strength of localization, i.e. the width of the plateaus 
closest to the CNP. Therefore, Eq. ([3]) cannot be used 
to relate the slope to the number of charge carriers. In 
fact, Eq. ([2]), which predicts a 1/n-behavior of the slope 
around the CNP when pB oo, anticipates such an 
overshoot. 

A similarly unexpected behavior for high temperature 
is observed in the resistance at high magnetic fields. 
Fig. mjc) shows the resistivity of BLG in the quantized 
regime as a function of the total charge carrier concen- 
tration for various temperatures. In contrast to the non- 
quantized regime, the resistivity exhibits a metal-like in- 
crease with temperature. 

As we will show in the following subsection, this pecu- 
liar transport properties in the quantized regime can be 
described in a simple DOS model with localized and ex- 
tended states where the temperature dependence of the 
resistivity is governed by a redistribution of charge car- 
riers due to the thermal activation. 



B. Model Calculations and Discussion 

The model used for our calculation starts with the 
quantized DOS for graphene in magnetic field consists 
of separate Gaussian-shaped Landau levels (LLs) as de- 
picted schematically in Fig. [SJ States in the middle of 
each LL are extended, and states are localized in the tails. 
For our simulations of the DOS in BLG shown below, we 
have used a constant width A = 100 K for the higher 
Landau levels and A — 200 K for the zero-energy LLJ^ 
Landau levels above (below) the CNP, populated by elec- 
trons (holes), have a degeneracy of AeB/h. The zeroth 
LL with a degeneracy of AeB/h for SLG and SeB/h for 
BLG is populated by electrons and holes simultaneously, 
such that electron and hole conduction can be found both 
above and below the CNP. We can then calculate the con- 
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FIG. 4: (color online) Transport coefficients pxy and pxx and of the BLG sample. Panel (a) and (c) show the experimentally 
measured values which are compared to the calculated quantities in (b) and (d). The insets in (a) and (c) show the corresponding 
slopes of the Hall resistivity at the CNP as a function of temperature. The insets in (b) and (d) show the temperature dependance 
of the conductivity at the CNP. 



ductivity, using the Kubo-Greenwood formalisnifiSi^: 

C^xx = CF^x + ^xx (7) 

dE. 



dE 



dE 



Here f{E) is the Fermi distribution and K'^-'^{E) is an 
energy dependent function containing the density of the 
extended states, the particle velocity and the scattering 
time, the upper indices e, h correspond to electrons and 
holes, respectively. Modeling K{E) by a superposition of 
Gaussians with the above mentioned widths, the Kubo- 
Greenwood formula reproduces qualitatively the behav- 
ior of the experimental conductivity in graphene, but 
does not give a universal value of minimal conductivity;^ 
Therefore, we normalized the integral in Eq. ([7]) such 
that, at zero temperature, axx at the GNP has the uni- 
versal value of Ae^/h in accordance with the experimen- 
tally observed value. 

Similarly, we can determine the Hall conductivity by 
summing up all extended states below the Fermi level, 
smeared out by temperature: 



/oo 
[D^,x,{l-f{E))-Dl,J{E)] dE 
-OO 



(8) 



When introducing a localized DOS, the calculated con- 
ductivity (T^y develops plateaus which are not at multi- 
ples of Ae^/h^ because the number of extended states 
within a LL is less than the total level degeneracy. In or- 
der to reproduce the correct number of occupied LLs in 
Eq. we have normalized the extended DOS for each 
LL to the total number of extended states within it. 

After calculating the conductivities we convert them 
into the resistivity pxx and the Hall resistivity pxy using 
the standard matrix relations, the results are shown in 
Fig. m^b) and (d). They can be compared to the exper- 
imental curves, plotted as panels (a) and (c) of the fig- 
ures. The calculated curves reproduce qualitatively the 
experimental results, and display a similar temperature 
dependence. In particular, at the CNP, the resistivity 
(insets in Fig. \M,c) and (d)) and the slope of the Hall 
resistivity (insets in Fig. |3](a) and (b)) increase with in- 
creasing temperature up to 200 K. Moreover, as is also 
seen in the experiment, the overshoots develop at the be- 
ginning of the plateaus at high temperatures (compare 
panels (a) and (b) in Fig. S]). 

Starting from about 500 K (as calculated for 30 T) 
the system smoothly enters the classical regime. The 
conductivity (inset to Fig. IH^d)) starts to increase with 
increasing temperature and the slope of the Hall resis- 
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FIG. 5: (color online) Schematic density of states (DOS) of 
broadened Landau levels in BLG. The dashed regions corre- 
spond to localized states between two levels. 

tivity (inset to Fig.|4l[b)) decreases. Experimentally, this 
trend in the non-quantized regime is verified by the low 
field data presented earlier in Figs. [2] and [3l 

V. CONCLUSIONS 

To conclude, we have investigated the coexistence of 
electron and hole magneto-transport in graphene for a 
wide range of temperatures and magnetic fields. For 



all temperatures and fields, the Hall resistivity smoothly 
crosses zero at the CNF indicating the simultaneous pres- 
ence of both electrons and holes at the CNF. In the 
non-quantized regime, the slope of the Hall resistivity 
decreases with increasing temperature, which is accom- 
panied by an increasing conductivity at the CNF. This 
behavior can be quantitatively modeled by a thermal ac- 
tivation of charge carriers experiencing relatively large 
potential fiuctuation of the order of 20 meV. 

In the quantum Hall regime, increasing the tempera- 
ture has the opposite effect on the slope of the Hall ef- 
fect and produces counter-intuitive overshoots when ap- 
proaching the v = ±2 plateaus. Using an appropriate 
density of states, we have shown in model calculations 
that such a behavior is expected for any two-carrier sys- 
tem with a quantized energy spectrum. 
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